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Abstract. The Mott metal-insulator transition in the two-band Hubbard model in infinite dimensions is 
studied by using the linearized dynamical mean-field theory recently developed by Bulla and Potthoff. 
The phase boundary of the metal-insulator transition is obtained analytically as a function of the on- 
site Coulomb interaction at the d-orbital, the charge-transfer energy between the d- and p-orbitals and 
the hopping integrals between p — d, d — d and p — p orbitals. The result is in good agreement with the 
numerical results obtained from the exact diagonalization method. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 71.27.+a Strongly correlated electron 
systems; heavy fermions - 71.30.+h Metal-insulator transitions and other electronic transitions 



1 Introduction 

The Mott metal-insulator transition (MIT) is a fundamen- 
tal problem in the theory of strongly correlated electrons. 
Recently, significant progress has been achieved in under- 
standing this transition by using dynamical mean-field 
theory (DMFT) §. In the DMFT, the lattice problem 
is mapped into an impurity problem embedded in an ef- 
fective medium by neglecting the momentum dependence 
of the self-energy. Various methods, such as the iterated 
erturbation theory jy, the non-crossing approximation 
the quantum Monte Carlo (QMC1 method ||, the ex- 
act diagonalization (ED) method [Q and the numerical 
renormalization group (NRG) method §,§, enable one 
to solve the corresponding impurity problem. The DMFT 
becomes exact in the limit of infinite spatial dimensions 
d = oo and is believed to be a good approximation in 
high dimensions. 

The Mott MIT in the half-filled single-band Hubbard 
model on the d — oo Bethe lattice is found to occur as 
a first-order phase transition below a critical temperature 
T c w 0.02W where W is the bare bandwidth jy. Below 
T c , a coexistence of the metallic and insulating solutions 
was found for the same value of the on-site Coulomb in- 
teraction U in the range U cl (T) < U < U c2 (T) (l],|,|,|l0], 
PUp!^ . At zero temperature, coexistence is also obtained 
for values of U such that U c \ < U < U C 2- The results of 
the ED method for the critical values of U (at T = ) ar e 
U cl w 1.2W and U c2 « 1.5W @ (see also Section g§. 
It agrees well with the recent NRG result U c \ s» 1.25VF 
and U C 2 ~ 1A7W (u|. The energy of the metallic state 
is lower than that of the insulating state for values of U 



in the range U cl < U < U c2 . Therefore the Mott MIT 
occurs at U — U C 2 as a continuous transition at T = 0. In 
this paper we will concentrate solely on the Mott MIT at 
T = and, so from this point onwards we will denote the 
critical value U C 2 simply by U c . 

The Mott MIT is observed in various 3d transition- 
metal compounds, which are classified into two types: the 
Mott-Hubbard type and the charge-transfer type |p^ , ^5| . 
In the Mott-Hubbard type such as Ti and V compounds, 
the d—d Coulomb interaction U is smaller than the charge- 
transfer energy A between d- and anion p-orbitals. In this 
case, the energy gap of the insulator is given roughly by 
U and a MIT occurs at a specific value of U as this inter- 
action strength is varied. In the charge-transfer type such 
as Co, Ni and Cu compounds, U is larger than A. Then 
the energy gap is roughly given by A and a MIT occurs at 
a critical value of A when this energy difference is varied. 

In the single-band Hubbard model, there is only the 
parameter U/W. The DMFT satisfactorily explains the 
Mott-Hubbard type MIT as a function of U/W. However 
we need at least a two-band Hubbard model with the pa- 
rameters U and A to describe the both types of MIT. 

In this paper we wish to study the Mott MIT with 
both the Mott-Hubbard type and the charge-transfer type 
mechanisms over the whole parameter regime. We use a 
two-band Hubbard model characterized by the following 
parameters: the on-site Coulomb interaction U at the d- 
orbital, the charge-transfer energy A between d- and p- 
orbitals and the hopping integrals t pc t, tdd and t pp between 
p — d, d — d and p — p orbitals, respectively. Several au- 
thors have studied the model using the DMFT approach 
j|,[l(],[lj],[l8, 19 1 . However, numerical problems make it dif- 
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ficult to obtain the critical values of the Mott MIT for this 
model, in contrast to the single-band Hubbard model. In 
the half-filled single-band Hubbard model on the Bethe 
lattice, the chemical potential is fixed to /i = ^ because 
of the particle-hole symmetry. On the other hand, for the 
two-band Hubbard model, the Mott MIT occurs away 
from particle-hole symmetry and so the chemical potential 
has to be determined explicitly to fix the electron density 
per unit cell to be unity. This calculation consumes a lot of 
CPU time. Furthermore, the Mott MIT point is a function 
of several parameters, all of which have to be calculated, 
making it a difficult numerical problem. 

In this work we show that there is an alternative ap- 
proach which has clear advantages as it can be handled 
analytically. It is based on the linearized version of DMFT, 
as developed by Bulla |2(J and applied by Bulla and Pot- 
thoff [|f]], where the critical value U c of the Mott MIT 
can be calculated analytically. In the linearized DMFT, 
the hybridization function between the impurity level and 
the conduction band is approximated by a single pole at 
the Fermi level. For the single-band Hubbard model on the 
Bethe lattice, the critical value is given by U c — 1.5W. The 
result is in good agreement with the numerical result from 
the ED and the NRG calculations of the full DMFT men- 
tioned above. The generalization of the linearized DMFT 
to more complicated lattices was also discussed in ref. . 
However, the models considered there were restricted to 
the particle-hole symmetric case where the chemical po- 
tential is fixed to the value fj, = 

In the present paper, we study a form of the two-band 
Hubbard model, which shows the Mott MIT away from 
the particle-hole symmetry. We generalize the linearized 
DMFT to the particle-hole asymmetric case and obtain 
an analytical expression for the critical values of the Mott 
MIT. A detailed account of the calculations is given in 
Section and the analytical results are compared with 
the numerical results from the ED method for the several 
values of parameters in Section 0. A good agreement be- 
tween the two approaches is found in every case. Some lim- 
iting cases in the Mott-Hubbard type and charge-transfer 
regimes are discussed in Section |^ and our conclusions are 
given in Section || 

2 Linearized dynamical mean-field theory 

2.1 Single-band Hubbard model 

First, we consider the single-band Hubbard model, 

H = - J2 l ^A c \a c jo + h.c.) + 4t c n c Ii c a- (!) 

<i,j>,(T i 

In the limit of infinite dimensions, the self-energy becomes 
purely site-diagonal and the DMFT becomes exact. The 
local Green's function G{z) can be given by the impurity 
Green's function of an effective single impurity Anderson 
model, 



+ E £fcC L^a +^14(/t Cfcff +4/,), (2) 

where Sf is the impurity level and £k are energies of con- 
duction electrons hybridized with the impurity by Vfc. In 
the model eq. (|2|), the non-interacting impurity Green's 
function, 

g {z) = (z-e f -A{z))-\ (3) 
with the hybridization function, 

^) = £-^L ,4) 

fe 

includes effects of the interaction at all the sites except 
the impurity site and is determined self-consistently so as 
to satisfy the self-consistency equation. 

For simplicity, the calculations in this paper are re- 
stricted to the Bethe lattice with the connectivity q and 
the hopping tjj = [p2| . In the limit q — oo, the self- 
consistent equation is given by 

Go(z)- 1 = z + ^-t 2 G(z), (5) 

where /x is the chemical potential for the original lattice 
model. In the non-interacting case, the local Green's func- 
tion is obtained from eq.(||) with Qo(z) = G(z). It yields 

a semicircular density of states: D(e) — ^ -y/f — (^j^) 2 

for \e + fi\ < 2t. Because of the particle-hole symmetry at 
half-filling, the chemical potential and the impurity level 
are set to fi = ^ and Ef — — ■j, respectively. Then, the 
self-consistency equation (||) is simply written by 

A(z) = t 2 G(z). (6) 

When the system approaches the MIT from the metal- 
lic side at T = 0, the central quasiparticle peak is found to 
appear to be isolated from the upper and the lower Hub- 
bard bands The width of the quasiparticle peak 
vanishes in the limit U — ► U c . In this limit, the effect 
of the Hubbard bands on the quasiparticle peak becomes 
rather small f^lf . Therefore, Bulla and Potthoff [^1) used 
an approximate form for the hybridization function where 
the contribution from the Hubbard bands are completely 
removed and that the quasiparticle peak is replaced by a 
single pole at the Fermi level p^| , 

with the small weight Aq which will be determined self- 
consistently. This model with eq.(0) corresponds to the 
two-site Anderson model , 

+£c Y J cW+vY,uio rT +ciui (8) 

(7 (7 
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with V = y/Ao, e c — and £/ — — The model is solved 
analytically to obtain the impurity Green's function which 
has four poles: two poles at ui rs ±¥ and two poles near 
the Fermi level ui « when U is large. These latter poles 
are precursors of the Kondo resonance in the Anderson 
model with complete conduction band. 

When U — > U c , the hybridization becomes V = y/ Aq — > 
0. In this limit, the impurity, and therefore, the local 
Green's function is given by G(z) — w/z near the Fermi 
level with the residue |24|] (see also Appendix |A|) , 



(9) 



up to the second order in V. From eqs. (0) and (||), we 
obtain a new hybridization function which has a pole at 
z = with the weight, 



A' 



(10) 



When we solve the self-consistency equation (g) by itera- 
tion, Aq for the (JV + l)th iteration step is expressed in 
terms of Aq for the TVth step through eq.([To|). Therefore 
the critical value for the MIT is given by 



U c = 6t = 1.5 W, 



(11) 



with the bare bandwidth W = At. For U < U c , the weight 
Aq increases exponentially with iteration number and, 
then, the single pole approximation for A{z) breaks down. 
For U > U c , Aq decreases exponentially to obtain the self- 
consistent value Aq — corresponding to the insulating 
solution. 



2.2 Comparison with numerical methods 

Here we estimate the reliability of the linearized DMFT by 
comparing the analytical result of the critical value eq. ( |lll ) 
with the available numerical results. The effective single 
impurity Anderson model eq.((^) is approximately solved 
by using the exact diagonalization of a cluster model with 
finite system size n s (the ED method). The Wilson's 
numerical renormalization group method is also used to 
solve the model eq.(||) in the thermodynamic limit |||| 
(the NRG method). 

Fig. |l| shows the U dependence of the quasiparticle 
weight defined by Z = (1 - ^§^U=o) _ \ with the local 
self-energy £{z) = Qq(z)~ 1 — G(z) _1 , calculated from the 
NRG method || and the ED method. The NRG and 
the ED results agree very well over the whole range of 
[/-values. 

When we approach the MIT point from the metallic 
side, Z continuously becomes zero at a critical value U c 
as shown in Fig. [I]. Recently, Bulla |l3| obtained the pre- 
cise result of the critical value U c = 5.88i by using the 
NRG method. In the ED method, an extrapolation of the 
systems with up to n s — 11 yields the n s — > oo extrapo- 
lated value U c = 5.87t [^5|. The result from the linearized 



DMFT, cq.(prl), is in very good agreement with the NRG 
and the ED results. It also agrees well with the value 
of U c = 5.84i from the projective self-consistent method 
(PSCM) (| and with the value of U c = 6.04i obtained 
in the NRG calculations of Shimizu and Sakai [^6|. The 
iterated perturbation method, where the effective single 
impurity Anderson model eq.(0) is solved within the sec- 
ond order perturbation with respect to U, gives a larger 
critical value U c = 6.6i as compared to the other non- 
perturbative approaches. The random dispersion approx- 
imation (RDA) [^7j predicts a considerably lower critical 
value U c = 4.0<. The origin of this discrepancy is presently 
not clear. 



2.3 Two-band Hubbard model 

Next, we consider the two-band Hubbard model on a Bethe 
lattice with connectivity q, 

H =HjT. (dl Pj(7 +h.c.) + Uj24^iAldil 

V <i,j>,cr i 

+ — { d \A'° + h - c -) + £ dYd\ a d ia 

^ <i,i / >,<T i,cr 

+ — E (p]aPj'-r + h.C.)+e p J2p]aPj-r, (12) 
^ <j-'j'>,v j,cr 

where d\ a and pt are creation operators for an electron 
with spin a in the c£-orbital at site i and in the p-orbital 
at site j, respectively. t p d, t&d an d tpp are the hopping 
integrals between the nearest neighbour p — d, p — p and 
d— d orbitals, respectively. The charge-transfer energy A is 
defined by A = e p — Ed > [p8[ . In eq.(^2|), we assume that 
p- and d-orbitals are on different sub-lattices of a bipartite 
lattice, more explicitly, a Bethe lattice with connectivity 
q. In the limit q — oo, the self-consistency equations for 
the local Green's functions are given by p|,^6| 

Gq(z)- 1 = z + »-e d - t 2 pd G p (z) - t 2 dd G d (z), (13) 
Gp(z)- 1 ^z + fi-Ep- t 2 pd G d {z) - t 2 pp G p {z), (14) 

where G p (z) is the local Green's function for the p-electron 
and Gd(z) is that for the <i-electron which can be given by 
the impurity Green's function of an effective single impu- 
rity Anderson model eq.(|). 

The non-interacting local Green's functions are eas- 
ily obtained by solving the self-consistency equations jl^ , 
|l4|) with Qq{z) — Gd(z). They yield the densities of states 
which consist of two bands: a d-band with the bandwidth 
Wd ~ 4:t dd around the d-level and a p-band with the band- 
width W p w 4i pp around the p-level when t p d is small. 
When the electron density per unit cell n — 1, the ri- 
band is half-filled and the system is metallic as shown in 
Fig.g(a). 

In the presence of U, the system is found to show the 
Mott MIT @0|ijyij. When U < A (Mott-Hubbard 
type), the energy gap of the insulating state is approxi- 
mately given by U (Fig. ||(b)). On the other hand, when 
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Fig. 1. The quasiparticle weight Z in the single-band Hubbard model on the Bethe lattice as a function of U: a comparison 



between the NRG result [13 
and that from the PSCM 



and the ED result (the system size n s = 8). The critical value from the linearized DMFT, C/ C LDMF , 
are also shown. We set t = 1. 



(a) d-band metal (U=0) 



(b) Mott-Hubbard type insulator 



(c) charge-transfer type insulator 
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Fig. 2. Schematic figures of the density of states in the two-band Hubbard model: (a) the d-band metal with (7 = 0, (b) the 
Mott-Hubbard type insulator with U < A and (c) the charge-transfer type insulator with U > A [eSI. Ef is the Fermi level. 



U > A (charge-transfer type), it is approximately given 
by A (Fig. |2j(c) ) . When the system approaches the MIT 
from the metallic side, the central quasiparticle peak is 
found to be largely isolated from the upper and lower 
Hubbard bands in the Mott-Hubbard type as observed 
in the single-band Hubbard model. It is also found to be 
largely isolated from the p-band and the lower Hubbard 
band in the charge-transfer type jl7],[l8j. The width of the 
quasiparticle peak vanishes in the limit of the MIT point 
in both types. 

Now, we study the MIT of this system by using the 
linearized DMFT. We assume that in the limit of the MIT 
point the effect of both the Hubbard bands and the p- 
band on the quasiparticle peak is negligible. Then we use 
an approximate form for the hybridization function eq. (^) 
as used in the single-band Hubbard model. In this case, 
the effective single impurity Anderson model eq.(^|) with 
eq.(0) corresponds to the two-site Anderson model cq.(||) 
with V = \/Ao, e c = and £/ = Ed = £d — The 



local d-Green's function in this model is obtained by (see 
Appendix 



Gd(z) = h 



Z Z — £i Z — £4 

where the residue Wd is (up to the second order in V) 



(15) 



2e 2 d 2(e d + Uf e d {e d + U) 



, (16) 



and w\, W4, E\ and £4 are given in Appendix 

To calculate the local p-Green's function, we assume 
an approximate form [f29| 



G p (z) = ^ + 



w v 1 — w D 



(17) 
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with Sp = Sp — fi. By using eq.fllTp in eq. (|14|) together with 
cq.(15), we obtain w p up to the second order in V: 



t2 p d w d 



(- t2 pd 



(18) 



Substituting eqs.([lj| |T^) into eq.(|l3|) yields a new hy- 
bridization function which has a pole at z — with the 
weight A' = FAq with 

F(t p d, t P p, tdd, U, Ed, £p) 

5 5 4 



2e 2 



2(e d + Uf e d (e d + U) 



t 4 

L pd 



id 



Wd+U) 



+ tld 



(19) 



- 1 2 
pp 



Following the same argument discussed in Section ETL we 
have an equation to determine the MIT point within the 
linearized DMFT: 



(20) 



In eq. ( pCf ) , F includes the chemical potential y, (through 
e~d and SpJTln the two-band Hubbard model, the MIT oc- 
curs away from the particle-hole symmetry as shown in 
Fig. |[ Then we have to determine the chemical poten- 
tial explicitly to obtain the critical values of the MIT. In 
general, the chemical potential is determined so as to fix 
the electron density. In the linearized DMFT, we focus on 
the low-energy part of the Green's function and determine 
it self-consistently. On the other hand, the high-energy 
part, whose details are neglected, is not determined self- 
consistently. Then we fail to obtain the precise expression 
of the electron density calculated from the Green's func- 
tion within the linearized DMFT. However, in the next 
paragraph we show that we can use an alternative condi- 
tion to determine ll, based on the fact that at the MIT 
point A' has a minimum value as a function of ll. This 
condition gives 



, tpp, tdd: i &di £p) 0* 



(21) 



Combined use of eq. (poj) and eq.(^l|) with eq.([l9|) enables 
us to obtain an analytic expression for the phase boundary 
separating the metallic and insulating regimes as a func- 
tion of the parameters: t p d, t pp , tdd, U and A = e v — Ed- 

In the metallic regime, the chemical potential Li(n) is 
continuous at n = 1 as a function of n. On the other 
hand, in the insulating regime, fJ.(n) has a jump at n = 1. 
When we approach the MIT phase boundary from the 
metallic side, n(n) is still continuous even in the limit of 
the MIT point. Correspondingly, there are three cases in 
the ll dependence of A (fi) as below. (1) In the metallic 
regime, A (fi) > A for all ll resulting in the metallic 
solution for all n. (2) In the insulating regime, A' (ll) < A 
for /i_ < /i < while, A' (ll) > Aq for ll < //_ or 



ll > Then the system is a Mott insulator for /i_ < 
ll < and /i shows a jump from /x_ to at n = 1. 
(3) On the phase boundary of the MIT, A' (p) = A for 
ll = Li(n = 1), while, A' (li) > Aq for ll ^ (j,(n = 1). 
Then A' q (ll) has a minimum at fi = fi(n = 1). Therefore 
the equation (|2l]) is the unique condition to determine the 
chemical potential on the MIT phase boundary within the 
linearized DMFT. 

In the single-band Hubbard model on the Bethe lattice, 
the condition to minimize A' is written by Qjjf- = where 
Wd is defined by eq.([l6|) with Ed — (Ed — — /i) without 
assuming the particle-hole symmetry. It yields the exact 
value of the chemical potential, /i = S, as expected. In 
the two-band Hubbard model, the chemical potential thus 
obtained agrees well with that from the ED method where 
the high-energy part is also determined self-consistently. 
It will be shown in Section [| 



3 Discussion 



3.1 Phase diagram 



From the coupled equations fl20|j2l| ) with eq.(|l9|), we can 
easily obtain the phase boundary of the MIT as a function 
of tpd, tpp, tdd, U and A within the linearized DMFT. 



Figs. |3|(a)-(c) show the phase diagrams of the MIT on the 
A — U plane for several values of t pp and tdd with t p d — 1 
together with the result obtained from the ED method for 
The result from the linearized DMFT 



tpp — tdd — [19 



is in good agreement with the ED result for all values of 
A and U in the case with t pp = tdd = 0. 

The MIT is observed when U is varied for U < A 
(Mott-Hubbard type), while it is observed when A is var- 
ied for U > A (charge-transfer type) as seen in Figs. ||(a)- 
(c). The phase boundary smoothly connects the Mott- 
Hubbard type and the charge-transfer type transitions for 
all values of tdd and t pp . As tdd and/or t pp increase, the 
metallic region monotonically increases. However, the ef- 
fect of tpp on the critical value of the Mott-Hubbard type 
transition is rather small as seen in Fig. ||(b) . This will be 
discussed in Section 

In the presence of t pp and tdd, we also calculated the 
quasiparticle weight Z by using the ED method. When we 
approach the MIT point from the metallic side, Z contin- 
uously becomes zero as shown in the insets of Figs. |i](a) 
and (b). The critical value U c (A c ) thus obtained is plot- 
ted in Fig. |(a) (Fig. |(b)) as a function of tdd (tpp) for 
a fixed value of A (U) together with that obtained from 
the linearized DMFT. The agreement between the two 
methods is good even in the case with finite t pp and tdd- 
Thus we conclude that the linearized DMFT gives a reli- 
able estimate for the phase boundary of the Mott MIT in 
the two-band Hubbard model over the whole parameter 
regime. We note that the chemical potential at the MIT 
point from the linearized DMFT also agrees well with that 
from the ED method. This confirms that eq.(pl|) is a reli- 
able condition to determine the chemical potential at the 
MIT point. 
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Fig. 3. Phase diagrams of the two-band Hubbard model at T = and n = 1. Solid lines are phase boundaries separating 
the metallic and insulating regimes obtained from the linearized DMFT as functions of A and U for several values of tdd with 
t pp = (a), for several values of t pp with tdd = (b) and for several values of t pp = 5tdd (c). Closed circles are the critical values 
for tpp = tdd = calculated from the exact diagonalization method ]l9(] . We set t p d — 1 in all figures. 



By eliminating the chemical potential /i directly from 
the coupled equations (^,^l|) with cq . (|l9|) , an analytic ex- 
pression for the phase boundary of the MIT is obtained 
within the linearized DMFT, although it is rather compli- 
cated. We can however get some simple analytical expres- 
sion in limiting cases as below. 



3.2 Mott-Hubbard regime (U < A) 

For the case with U < A, the MIT is observed when U 
is varied as seen in Figs. g(a)-(c). In the limit A —> oo, 
eqs.( [20|pl| ) with cq.(|l9|) yield the chemical potential [i w 
and the critical value U c : 



i 4 

6\ -Ei +t 2 



When t pd = 0, eq.fl22|) results in U c ~ 6t dd — 1.5Wd with 
the bare d-bandwidth Wd — 4£dd which is equivalent to 
cq.(|TT|) in the case with sin gle- band Hubbard model as 
expected. When t dd = 0, eq.p2]) yields U c ~ 1.5W d with 
the bare d-bandwidth of the d-p hybridized band given 



4t„ 



by W d = \[{A 2 + 16*2 ) 3 - A] « In general, U c is 
approximately given by 1.5 times the bare c?-bandwidth 
in the Mott-Hubbard regime. We note that, in the limit 
A — > oo, U c is independent of t pp as shown in eq.(|22]) and 
also shown in Fig. ^(b). This enables us to describe the 
MIT in the Mott-Hubbard type by using the single-band 
Hubbard model. 



3.3 Charge-transfer regime (U > A) 



For the case with U > A, the MIT is observed when A 
' is varied as seen in Figs. @(a)-(c). In the limit U — » oo, 
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L dd l pp 

Fig. 4. The critical value U c as a function of tdd for t pp = and A = 10 (a), and the critical value A c as a function of t pp for 
tdd = and (7 = 6 (b), obtained from the linearized DMFT (solid line) and from the ED method (closed circle). Insets show 
the quasiparticle weight Z from the ED method (the system size n a = 8) as a function of U for several tdd with t pp = and 
A — 10 (a), and as a function of A for several t pp with tdd = and U — 6 (b). We set t p d = 1 in all figures. 



eqs. (|2Qj2l| ) with eq. © yield the critical value A c 



4 Conclusions 



2(\/i-! 

l + \/5£ 



2_ "PP 



tpd -*^ > tppi^ddi 
tpp tpdftddt 
tdd tpd-, tpp. 



(23) 



Correspondingly, the chemical potential is given by /i ss 
£p ^ £d for t pti > t P p, tdd, V ~ £d for t pp > t pd ,t dd and 
/i w e p for tdd S> tpd, tpp- We note that for the small value 
of A < A c the system is metallic even in the limit U — > oo 
as mentioned in ref. jlB). When U » A ^> t p d,t pp ,tdd, 
d-orbitals are almost singly occupied and p-orbitals are 
nearly empty. The electron transfer from a <i-orbital to 
a p-orbital costs the charge-transfer energy A while it 



gains the kinetic energy: K 



tpd s^-" tpp, tdd-, 



K 



tpp for tpp 



> t pd Am and K - t dd + O(-ff-) for 
tdd ^> tpd, tpp. Equation (|23| ) means an energy balance be- 
tween the charge-transfer energy and the kinetic energy 
& 

For general values of tpd,t pp and tdd, the explicit ex- 
pression for A c is still complicated even in the limit U — * 
oo. However the critical value A c is roughly given by the 
energy balance mentioned above in the charge-transfer 
regime. We note that A c is positive for all values of U, 
tpd, tpp and td d in contrast to the negative- A c predicted 
by the local impurity approximation |3l| where the im- 
purity Anderson model is solved (not self-consistently) to 
determine the critical value of the MIT. This descrepancy 
will be discussed in Section 13. 



Within the linearized DMFT, we have obtained a good 
description of the phase diagram of the Mott MIT analyt- 
ically in the two-band Hubbard model over the whole pa- 
rameter regime including the Mott-Hubbard regime, charge- 
transfer regime and the intermediate regime. The ana- 
lytical result agrees well with the numerical result ob- 
tained from the ED method. Although the ED method is 
an approximate calculation, it is in very good agreement 
with the recent NRG method in the single-band Hubbard 
model and is expected to be reliable for the two-band Hub- 
bard model. 

We have used the two-site Anderson model in the lin- 
earized DMFT. The same model was also used in the ED 
method with the smallest cluster size n s = 2 and obtained 
the self-consistent solution which does not show the MIT 
up to U = lOOt. The difference between two methods is 
the self-consistent procedure. In the linearized DMFT, we 
focus only on the lowest-energy poles and determine them 
sclf-consistently. On the other hand, in the ED method, 
all of the poles are determined self-consistently so as to 
satisfy the self-consistency condition as close as possible 
even in the case with n s = 2. The n s = 2 system is insuffi- 
cient to describe the high-energy poles and fails to obtain 
the MIT within the ED method, while it is sufficient to 
describe the low-energy poles to obtain the critical value 
of the MIT within the linearized DMFT. 

In the linearized DMFT, we have assumed that as the 
MIT is approached the central quasiparticle peak becomes 
isolated from the upper and lower Hubbard bands and 
the p-band. Such an isolated quasiparticle peak occurs 
for the single-band Hubbard model p|,^3[ and also for 
the two-band Hubbard model jl7],[l^]. The phase bound- 
ary smoothly connects the Mott-Hubbard type and the 
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charge-transfer type transitions, showing that the cen- 
tral quasiparticle peak also smoothly changes from one 
type to the other. Even in the case of large p- and/or d- 
bandwidth, a finite t p d yields a hybridization gap between 
the p- and the d-bands and a central quasiparticle peak 
is found within the hybridization gap isolated from these 
bands in the limit of the MIT point. 

When tpd = and U > U c ~ 6idd, there is a transition 
from the Mott insulator to a metal at A c = 2t pp = -s 2 -. 
In this case the p and the d electron states are decoupled. 
There is a correlated d-band and a separate free p-band, 
and the transition at A c — 2t pp is simply due to the over- 
lap of these bands. The quasiparticle weight in the p-band 
remains unity as the transition is approached. This is quite 
different from the Mott MIT considered here, which is a 
many-body transition where the weight and the width of 
the quasiparticle peak within the hybridization gap due to 
tpd decrease to zero. The value of A c for this transition, 
given by eq.(|23|) in the limit t p d — > 0, differs from the 
value A c 
that tpd 
model. 

Zaanen, Sawatzky and Allen previously obtained a sim- 
ilar phase diagram of the Mott MIT by using the local im- 
purity approximation (LIA) [ jl5| . They mapped the lattice 
model onto the impurity Anderson model where the hy- 
bridization function is assumed and is not determined self- 
consistently in contrast to the DMFT. In the DMFT, the 
self-consistency condition plays an important role to take 
into account of the translational symmetry. By neglect- 
ing the translational symmetry, the LIA fails to describe 
precisely the metallic state and the MIT phase boundary 
although it gives reasonable description for the insulating 
state when the gap is large j32|. For the case with small 
p-bandwidth W p = 4i pp , the LIA predicts a negative value 
of A c in the charge-transfer type J3l| in contrast to the 
DMFT where A c is always positive. As mentioned in Sec- 
A c is roughly given by an energy balance between 



= 2t pp for the simple overlap transition, showing 
is a singular point of the two-band Hubbard 



tion 3.3 



the charge-transfer energy and the kinetic energy. Even in 
the case with t pp — tdd = 0, the kinetic energy due to t p d, 
which is underestimated in the LIA, gives a positive value 
of A c in the DMFT. 

Within the linearized DMFT, we have determined the 
chemical potential by using the condition that the hy- 
bridization function has a minimum instead of fixing the 
electron density to be unity. The chemical potential thus 
obtained is equivalent to the exact result in the single- 
band Hubbard model on the Bethe lattice. It also agrees 
well with the result from the ED method in the two- 
band Hubbard model. This method is applicable for the 
single-band Hubbard model without particle-hole symme- 
try such as the fec-type lattice in d = oo jj3) . It is also in- 
teresting to apply this method to more complicated mod- 
els, e.g., the model with the orbital degeneracy and the 
Hund rule coupling, where the numerical method of the 
full DMFT becomes rather difficult. 
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Appendix A 

Here we discuss the two-site Anderson model eq.(|^) p4| . 
We assume that the conduction level is between the atomic 
/-level and the upper Hubbard level: Ef<e c <£f + U. 
The one electron eigenstates 



\E±) = a±/+|0)+#fcC+|Q>, 
correspond to the eigenenergies 



(24) 



I 



E±= g {e c + £f±^(£c-£f) 2 +W2) . (25) 

For the small hybridization V 2 -C e c — £f , eq. ( p5| ) is sim- 
plified as 

V 2 V 2 
E+=e c + , £-=£/ , (26) 

E c ~ £f E c - Ef 

to leading order in V 2 , with the corresponding eigenstates 

V 



E C - Ef 



E,) = a 



\E-)=alf+- 



V 



Ec-Ef 



a ] v J I 3 



|o>, 



(27) 



(28) 



v 

with a = 1 — tjt — - — Tj- . Similarly, we obtain the three 
electron (one hole) eigenenergies 



E± = i Ue c + 3e/ + U± ^(£f + U-£ c ) 2 +4V 2 j(29) 

For the small hybridization V 2 <C £/ + U — e c , cq.(]2^) is 
simplified as 



E+ =£ c + 2£ f + U 
E_ = 2£ c + £f 



V 2 



Ef + U - Ec 

V 2 



Ef + U - Ec 



(30) 
(31) 



to leading order in V 2 . The corresponding eigenstates are 



\E+) = a 



V 



Ef + U-Ec 



U + c a |4>, (32) 



\E-) = a[U 



V 

Ef + U - £ c 



Ca |4), (33) 
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with a = 1 - 2(ef+ ^_ Ee)a and |4) = /+/+C+C+ 10). 

The two electron states can be classified as singlets 
or triplets. In the triplet state, the spatial part of the 
wavefunction is antisymmetric and the interaction U plays 
no role. Then the total energy of the triplet state is given 
by E+ + E- = e c + £/■ There are three possible singlet 
states which can be written by the linear combination of 
the states, 



1 



|<Ai) = - 7 =(cf/+-c+/f)|0) 

|02> =C+C+|0), 

|<A 3 ) = / T + / i + |0). 



(34) 

(35) 
(36) 



The eigenenergies are given by the solutions of the equa- 
tion, 



E-E c -ef -V2V -V2V 
-V2V E-2e c 
-V2V E-2e f -U 



= 0. 



To leading order in V 2 , the eigenenergies are 



E 1 = £ c + £f 

E 2 = 2e c + - 



2V 2 



2V 2 



£c 
2V 2 



E 3 = 2 £f + U + 



£f £f + U 



2V 2 



£f + U - £ c 

and the corresponding eigenstates are 




with a\ = 1- 



v 2 - v 



(37) 

(38) 
(39) 
(40) 



I0 3 ) (,41) 



(42) 

(43) 
and 



i ; = 1 — ( g/+ ^_ £ )2 ■ In our situation with Sf < £ c < £f + 



U, we find the singlet ground state, \Ei), with an energy 



gain of 



2V A 



2V Z 



due to the hybridization. This is 



-Ef Sf+U— £ e 

equivalent to the energy gain 2 J in the s-d model with the 
Kondo coupling J corresponding to the Schricffer- Wolff 
transformation for V 2 <C £ c — £f and V 2 <C £/ + U — e c . 

Now we calculate the /-electron Green's function of 
this model. When a /, ] electron is removed from the 
ground state \E\), there are two possible final states: \E + ) 
and \EJ). Correspondingly, there are two possible single- 
hole excitations with excitation energies, 



3V 2 

E + -E x = -£ f + + 



2V 2 



E— — E\ 



£ c ~£f £f + U-£ c 

V 2 2V 2 
-£c + + 



£ c — £f £f + U — £< 



= -£i, (44) 
= ~£2, (45) 



to leading order in V 2 . The matrix elements for these tran- 
sitions are 



(E + \f,\ El ) = ^[l 



2V 2 



(E_\f^\E 1 ) = ^l 



V2 V {£c-£f){£f + U-£ c )J ' 
2V 



y/2 V £ c - £f £f + U - £ C J ' 
which yield the transition probabilities: 

1 3V 2 V 2 



l(£+l/ T |£i 



2 2{£ C - £f) 2 (£f +U- £ C ) 2 

2V 2 



{£ c ~£f)(£f + U-£ C ) 



= lOi, 



(46) 



l(i?-l/Tl^i)| 2 



+ 



2V 2 



2{£ C - £f) 2 (£f +U- £c) 2 

2^2 

{£ c -£f)(£f + U -£ c ) 



= W 2 , 



(47) 



to leading order in V 2 . 

When a /, | electron is added to the ground state \E\), 
possible final states are \E_) and \E + ). Correspondingly, 
there are two possible single-particle excitations with ex- 
citation energies, 



- E l = £ C + 



2V 2 



V 2 



E+-E 1 =£ f + U + 



£c - £f Ef + U — Ec 

2V 2 W 2 



£3, (48) 



Ec - £f Ef + U - Ec 

= £4, (49) 

to leading order in V 2 . The matrix elements for these tran- 
sitions are 



aa\ ( 2V 



+ 



V 



£ c -£/ £f + U-E c 
2V 2 



T '^ i/ V2 V (e c -£/)(£/ + U-s c )J' 
which yield the transition probabilities: 

2V^_ 

(e c -£/)(£/ + U -E C ) 

= w 3 , (50) 
1 V 2 "W 2 



2 (E c -E f ) 2 2(£ f + U-£c) 2 

2V^ 

(e c -£/)(£/ + U -£ c ) 



(51) 



to leading order in V 2 . 
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From eqs.(|44|)-(|5l]), we obtain the /-electron Green's 
function which has four poles: 



W; 



Z - £i 



(52) 



In the limit V — > 0, high-energy poles at Ej rs and 
£4 ~ Sf + U have large residues w% w W2 ~ §, while low- 
energy poles merge together at £2 ~ £3 ~ with small 
total residue to = 102 + 103: 



5^ 2 



+ 



5^ 2 



2(e c - e/ ) 2 2( e/ + C/- £c ) 2 
(e c -e/)(e/ + £/-£ c )' 



(53) 
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